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Abstract 



> . 

■ We consider an 0{N) model coupled self-consistently to gravity in the semi- 

. classical approximation, where the field is subject to 'new inflation' type initial 

0^ , conditions. We study the dynamics self-consistently and non-perturbatively 

with non-equilibrium field theory methods in the large N limit. We find that 
0^ I spinodal instabilities drive the growth of non-perturbatively large quantum 

fiuctuations which shut off the infiationary growth of the scale factor. We 
find that a very specific combination of these large fluctuations plus the in- 
Oh. flaton zero mode assemble into a new effective field. This new field behaves 

r~| , classically and it is the object which actually rolls down. We show how this 

reinterpretation saves the standard picture of how metric perturbations are 
generated during inflation and that the spinodal growth of fluctuations dom- 
inates the time dependence of the Bardeen variable for superhorizon modes 
during inflation. We compute the amplitude and index for the spectrum of 
scalar density and tensor perturbations and argue that in all models of this 
type the spinodal instabilities are responsible for a 'red' spectrum of primor- 
dial scalar density perturbations. A criterion for the validity of these models is 
provided and contact with the reconstruction program is established validat- 
ing some of the results within a non-perturbative framework. The decoherence 
aspects and the quantum to classical transition through inflation are studied 
in detail by following the full evolution of the density matrix and relating the 
classicality of cosmological perturbations to that of long-wavelength matter 
fluctuations. 



*Laboratoire Associe au CNRS UA280. 
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I. INTRODUCTION AND MOTIVATION 



Inflationary cosmology has come of age. From its beginnings as a solution to purely 
theoretical problems such as the horizon, flatness and monopole problems it has grown 
into the main contender for the source of primordial fluctuations giving rise to large scale 
structure There is evidence from the measurements of temperature anisotropics in the 

cosmic microwave background radiation (CMBR) that the scale invariant power spectrum 
predicted by generic inflationary models is at least consistent with observations [^|-|^ and we 
can expect further and more exacting tests of the inflationary power spectrum when the MAP 
and PLANCK missions are flown. In particular, if the fluctuations that are responsible for 
the temperature anisotropics of the CMB truly originate from quantum fluctuations during 
inflation, determinations of the spectrum of scalar and tensor perturbations will constrain 
inflationary models based on particle physics scenarios and probably will validate or rule 
out speciflc proposals |B|,[7|,ITD|,ITT| . Already current bounds on the spectrum of scalar density 



perturbations seem to rule out some versions of 'extended' inflation |T0 . 

The tasks for inflationary universe enthusiasts are then two- fold. First, models of inflation 
must be constructed that have a well-deflned rationale in terms of coming from a reasonable 
particle physics model. This is in contrast to the current situation where most, if not all 
acceptable inflationary models are ad-hoc in nature, with flelds and potentials put in for the 
sole purpose of generating an inflationary epoch. Second, and equally important, we must be 
sure that the quantum dynamics of inflation is well understood. This is extremely important, 
especially in light of the fact that it is exactly this quantum behavior that is supposed to give 
rise to the primordial metric perturbations which presumably have imprinted themselves in 
the CMBR. This latter problem is the focus of this paper. 

The inflaton must be treated as a non- equilibrium quantum fleld . The simplest way to 
see this comes from the requirement of having small enough metric perturbation amplitudes 
which in turn requires that the quartic self coupling A of the inflaton be extremely small, 
typically of order ~ 10~^^. Such a small coupling cannot establish local thermodynamic 
equilibrium (LTE) for all fleld modes; near a phase transition the long wavelength modes 
will respond too slowly to be able to enter LTE. In fact, the super horizon sized modes will 
be out of the region of causal contact and cannot thermalize. We see then that if we want 
to gain a deeper understanding of inflation, non-equilibrium tools must be developed. Such 
tools exist and have now been developed to the point that they can give quantitative answers 
to these questions in cosmology ||T^16|. These methods permit us to follow the dynamics 



of quantum flelds in situations where the energy density is non-perturbatively large (~ 1/A). 
That is, they allow the computation of the time evolution of non-stationary states and of 
non-thermal density matrices. 

Our approach is to apply non-equilibrium quantum fleld theory techniques to the situa- 
tion of a scalar fleld coupled to semiclassical gravity, where the source of the gravitational 
field is the expectation value of the stress energy tensor in the relevant, dynamically chang- 
ing, quantum state. In this way we can go beyond the standard analyses [|T^-|2^ which treat 
the background as fixed. 

We will mainly deal with 'new infiation' scenarios, where a scalar field (p evolves under 
the action of a typical symmetry breaking potential. The initial conditions will be taken so 
that the initial value of the order parameter (the field expectation value) is near the top of 
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the potential (the disordered state) with essentially zero time derivative. 

What we find is that the existence of spinodal instabilities, i.e. the fact that eventually 
(in an expanding universe) all modes will act as if they have a negative mass squared, drives 
the quantum fiuctuations to grow non-perturbatively large. We have the picture of an initial 
wave-function or density matrix peaked near the unstable state and then spreading until it 
samples the stable vacua. Since these vacua are non-perturbatively far from the initial state 
(typically ~ m/\/A, where m is the mass scale of the field and A the quartic self-coupling), 
the spinodal instabilities will persist until the quantum fiuctuations, as encoded in the equal 
time two-point function ($(a;, t)^), grow to 0(m^/A). 

This growth eventually shuts off the inflationary behavior of the scale factor as well as 
the growth of the quantum fluctuations (this last also happens in Minkowski spacetime [|T^ ). 

The scenario envisaged here is that of a quenched or supercooled phase transition where 
the order parameter is zero or very small. Therefore one is led to ask: 

a) What is rolling down?. 

b) Since the quantum fluctuations are non-perturbatively large ( ~ ^/^), will not they 
modify drastically the FRW dynamics?. 

c) How can one extract (small?) metric perturbations from non-perturbatively large fleld 
fluctuations? 

We address the questions a)-c) as well as other issues below. 



II. NON-EQUILIBRIUM QUANTUM FIELD THEORY, SEMICLASSICAL 

GRAVITY AND INFLATION 

Our program consists of flnding ways to incorporate the non-equilibrium behavior of the 
quantum flelds involved in inflation into a framework that treats gravity self consistently, 
at least in some approximation. We do this via the use of semiclassical gravity |2l|| where 



we say that the metric is classical, at least to flrst approximation, whose source is the 
expectation value of the stress energy tensor (T^^) where this expectation value is taken 
in the dynamically determined state described by the density matrix p{t). This dynamical 
problem can be described schematically as follows: 

1. The dynamics of the scale factor a{t) is driven by the semiclassical Einstein equations 
1 

-G^i, + Qf^y + (higher curvature) = -{T^u)r- (2.1) 



Here Gr, Ar are the renormalized values of Newton's constant and the cosmological 
constant, respectively and G^u is the Einstein tensor. The higher curvature terms 
must be included to absorb divergences. 

2. On the other hand, the density matrix p{t) that determines (T^i/)_R obeys the Liouville 
equation 

^^ = [H,pit)], (2.2) 
where H is the evolution Hamiltonian, which is dependent on the scale factor, a{t). 



3 



It is this set of equations we must try to solve; it is clear that initial conditions must be 
appended to these equations for us to be able to arrive at unique solutions to them. Let us 
discuss some aspects of the initial state of the field theory first. 



A. On the initial state: dynamics of phase transitions 

As we mentioned above, the situation we consider is one in which the theory admits a 
symmetry breaking potential and in which the field expectation value starts its evolution 
near the unstable point. There is an issue as to how the field got to have an expectation 
value near the unstable point (typically at $ = 0) as well as an issue concerning the initial 
state of the non-zero momentum modes. The issue of initial conditions is present in any 
formulation of inflation but chaotic. 

Since our background is an FRW spacetime, it is spatially homogeneous and we can 
choose our state p{t) to respect this symmetry. Starting from the full quantum field ^{x,t) 
we can extract a part that has a natural interpretation as the zero momentum, c-number 
part of the field by writing: 

<l>(f,t) = (j){t) + 

0(t)=Tr[p(t)$(f,t)] = ($(f,t))). (2.3) 

The quantity \l/(x, t) represents the quantum fluctuations about the zero mode (j){t) and 
clearly satisfies (\&(x, t)) = 0. 

We need to choose a basis to represent the density matrix. A natural choice consistent 
with the translational invariance of our quantum state is that given by the Fourier modes, 
in comoving momentum space, of the quantum fluctuations \E'(x, t): 

^exp(-z /c-f) Mt)- (2.4) 

In this language we can state our ansatz for the initial condition of the quantum state 
as follows. We take the zero mode (j)(t = 0) = (po, (p(t = 0) = 0, where 0o will typically 
be very near the origin, while the initial conditions on the the nonzero modes ipk{t = 0) 
will be chosen such that the initial density matrix p{t = 0) describes a vacuum state (i.e. 
an initial state in local thermal equilibrium at a temperature Tj = 0). There are some 



subtleties involved in this choice. First, as explained in |T6[, in order for the density matrix 
to commute with the initial Hamiltonian, we must choose the modes to be initially in the 
conformal adiabatic vacuum (these statements will be made more precise below). This choice 
has the added benefit of allowing for time independent renormalization counterterms to be 
used in renormalizing the theory. 

We are making the assumption of an initial vacuum state in order to be able to proceed 
with the calculation. It would be interesting to understand what forms of the density matrix 
can be used for other, perhaps more reasonable, initial conditions. 

The assumptions of an initial equilibrium vacuum state are essentially the same used by 



Linde |T^, Vilenkin |T8| , as well as by Guth and Pi |^ in their analyses of the quantum 
mechanics of inflation in a fixed de Sitter background. 
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As discussed in the introduction, if we start from such an initial state, spinodal insta- 
bihties will drive the growth of non-perturbatively large quantum fluctuations. In order to 
deal with these, we need to be able to perform calculations that take these large fluctuations 
into account. Although the quantitative features of the dynamics will depend on the initial 
state, the qualitative features associated with spinodal instabilities will be fairly robust for 
a wide choice of initial states that describes a phase transition with a spinodal region in field 
space. 



III. THE MODEL AND EQUATIONS OF MOTION 

Having recognized the non-perturbative dynamics of the long wavelength fluctuations, 
we need to study the dynamics within a non-perturbative framework. That is, a frame- 
work allowing calculations for non-perturbatively large energy densities. We require that 
such a framework be: i) renormalizable, ii) covariant energy conserving, iii) numerically 
implementable. There are very few schemes that fulfill all of these criteria: the Hartree and 
the large approximation []TRO,|Il|. Whereas the Hartree approximation is basically a 
Gaussian variational approximation p2| , p3| that in general cannot be consistently improved 
upon, the large A^ approximation can be consistently implemented beyond leading order 
2^ , |25[] and in our case it has the added bonus of providing many light fields (associated 



with Goldstone modes) that will permit the study of the effects of other fields which are 
lighter than the inflaton on the dynamics. Thus we will study the inflationary dynamics of 
a quenched phase transition within the framework of the large A^ limit of a scalar theory in 
the vector representation of 0{N). 

We assume that the universe is spatially flat with a metric given by 



- a^{t) dx^ 



(3.1) 



The matter action and Lagrangian density are given by 



Sri 



d X 



d^x a^{t) 



2 ^ ' 2 a?{t) 



- V{^{x)) 



(3.2) 



2A^m2' ^ 




v{^) = — \^' — +^7^<^)^ (3.3) 



-7?m a('<^^t) d\t)\ 

where we have included the coupling of $(x) to the scalar curvature 'R-{t) since it will arise 
as a consequence of renormalization [|1^. 

The gravitational sector includes the usual Einstein term in addition to a higher order 
curvature term and a cosmological constant term which are necessary to renormalize the 
theory. The action for the gravitational sector is therefore: 
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d'^x a^{t) 



(3.5) 



with K being the cosmological constant (we use K rather than the conventional A/8ttG to 
distinguish the cosmological constant from the ultraviolet cutoff A we introduce to regular- 
ize the theory; see section IV). In principle, we also need to include the terms R^'^R^^ and 
aPfiu as they are also terms of fourth order in derivatives of the metric (fourth adia- 
batic order), but the variations resulting from these terms turn out not to be independent 
of that of 7^2 in the flat FRW cosmolo gy we are considering. 

The variation of the action S = Sg + Sm with respect to the metric g^i, gives us Einstein's 
equation 



G 



8nG 



+ aHau + Kg 



-<T.,> , 



(3.6) 



where G^i, is the Einstein tensor given by the variation of y/—gTZ, H^i, is the higher order 
curvature term given by the variation of y/—gTZ'^, and T^^ is the contribution from the 
matter Lagrangian. With the metric ( p.l| ), the various components of the curvature tensors 
in terms of the scale factor are: 



^0 



-3(a/a)^ 



^ a a \ 

= -6 - + - ' 

-6 -7^ + — 7^ - —V} 



12 



-6 (7^ + 3-7^ 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



Eventually, when we have fully renormalized the theory, we will set = and keep as our 
only contribution to Kji a piece related to the matter fields which we shall incorporate into 
T, 



A. The Large N Approximation 

To obtain the proper large limit, the vector field is written as 

$(x, t) = ((t(x, t), 7r(x, t)), 
with 7? an (A^ — l)-plet, and we write 

a(f,t) = ViV0(t) + x(x,t) ; (a(f,t)) = v^0(t) ; (x(f,t)) = 0. 



(3.11) 



To implement the large A^ limit in a consistent manner, one may introduce an auxiliary field 
as in ||25|. However, the leading order contribution can be obtained equivalently by invoking 



the factorization |IT^,|T6ll: 
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X 



[IT • TC) 



-'2 2 
TT X 



6(x^)x^ + constant , 



(3.12) 
(3.13) 
(3.14) 
(3.15) 
(3.16) 



To obtain a large N limit, we define [0,0] 



7r{x,t) = V^(f,t)(l,l,---,1), 
where the large limit is implemented by the requirement that 

(^2)^0(1), (x')^O(l), 0^O(1). 



(3.17) 



(3.18) 



The leading contribution is obtained by neglecting the 0{1/N) terms in the formal limit. 
The resulting Lagrangian density is quadratic, with linear terms in x and vf. The equations of 
motion are obtained by imposing the tadpole conditions < x(a;, t) >= and < n{x, t) >= 
which in this case are tantamount to requiring that the linear terms in x ^ind vr in the 
Lagrangian density vanish. Since the action is quadratic, the quantum fields can be expanded 
in terms of creation and annihilation operators and mode functions that obey the Heisenberg 
equations of motion 



TT (x, t) 



(Pk 

(2^ 



aufkit) e^'-^ + dl flit) e 



-ik-x 



The tadpole condition leads to the following equations of motion 0,116 

4>{t) + 3H{t) <p{t)+M\t) (Pit) 

with the mode functions 



0, 



Mt) = 0, 



where 



M^t) = -m' + ^n + ^02(t) + ^ii^'it)) . 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



In this leading order in 1/A^ the theory becomes Gaussian, but with the self-consistency 
condition 



mt)) 



d'k |/,(t)P 
(27r)3 2 



(3.23) 



The initial conditions on the modes /^(t) must now be determined. At this stage it proves 
illuminating to pass to conformal time variables in terms of the conformally rescaled fields 
(see and section IX for a discussion) in which the mode functions obey an equation 



7 



which is very similar to that of harmonic oscillators with time dependent frequencies in 
Minkowski space-time. It has been realized that different initial conditions on the mode 
functions lead to different renormalization counterterms |T^; in particular imposing initial 
conditions in comoving time leads to counterterms that depend on these initial conditions. 
Thus we chose to impose initial conditions in conformal time in terms of the conformally 
rescaled mode functions leading to the following initial conditions in comoving time: 



fkito) 



fkito) 



a{to 



a{to 



(3.24) 



with 



W^ = e + M\to 



7^(to) 
6 



(3.25) 



For convenience, we have set ci(to) = 1 in eq. ( p.25| ). At this point we recognize that when 
-^^(^o) — T^{to)/Q < the above initial condition must be modified to avoid imaginary 
frequencies, which are the signal of instabilities for long wavelength modes. Thus we define 
the initial frequencies that determine the initial conditions ( 3.24| ) as 



wl = k'' + M\t 



2/ 



7^(^o) 



6 

7^(to) 

6 



for k < 



n{to) 



for > 



7^(^o 



6 



(3.26) 
(3.27) 



As an alternative we have also used initial conditions which smoothly interpolate from 
positive frequencies for the unstable modes to the adiabatic vacuum initial conditions defined 
by ( p.24| , p^25D for the high k modes. While the alternative choices of initial conditions 
result in small quantitative differences in the results (a few percent in quantities which 
depend strongly on these low-A; modes), all of the qualitative features we will examine are 
independent of this choice. 

In the large limit we find the energy density and pressure density to be given by P^JTB 



^ 1 ;2 
— = -0 H — m 

N 2^ 2 



o Za CL 



3p 



N 



-</)^ + 2m^ 0^ + - 



+ 



2m^ 



A 



a 

+ ^[202 (^2) + {^Y] + ((^ ^) + (^^) + 3 - m 

2 \ a 



a 
a 



(3.28) 



(3.29) 



where {ip"^) is given by equation ( p. 23 ) and we have defined the following integrals: 
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<Pk 



2(27r)3 



(3.30) 
(3.31) 



The composite operators (ipip) and (ipip) are symmetrized by removing a normal ordering 
constant to yield 



^m) + i^) = -J 



dt 



(3.32) 
(3.33) 



The last of these integrals, ( p.33| ), may be rewritten using the equation of motion (|3.21| ): 

(3.34) 



a a"^ 



It is straightforward to show that the bare energy is covariantly conserved by using the 
equations of motion for the zero mode and the mode functions. 



IV. RENORMALIZATION 



Renormalization is a very subtle but important issue in gravitational backgrounds ^1 
The fluctuation contribution (?/;^(a;, t)), the energy, and the pressure all need to be renor- 
malized. The renormalization aspects in curved space times have been discussed at length 
in the literature and have been extended to the large N self-consistent approximations 
for the non-equilibrium backreaction problem in p5|JT6|j26[] . More recently a consistent and 
covariant regularization scheme that can be implemented numerically has been provided 



In terms of the effective mass term for the large N limit given by (|3.22|) and defining the 
quantity 



B{t) = a\t) (M\t) - n/Q) , 



M'{t) = -ml + ^snt) + ^cp\t) + ^{i^\t))s 



(4.1) 
(4.2) 



where the subscript B stands for bare quantities, we find the following large k behavior for 
the case of an arbitrary scale factor a{t) (with a(0) = 1): 



1 



1 



+ 



ka'^it) 2k^a^{t) 
1 



Bit) 



d 



5(2) + 0(1/^) , 



Oil/k' 



(4.3) 



lA(t)p 



+ 



a^{t) 2ka\t) 



B{t) + 2d 
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+ 



5(1) + 0(1/^5 



{-B{ty - a{tfB{t) + ?,a(t)a(t)B{t) - Aa^{t)B(t)] + 



(4.4) 



fk{t)m+h{t)fm 



a{t) 



1 



ka'^{t)a{t) Ak^a'^it) 



ait) , , 
2-^ B{t) 
a{t) 



+ 0{l/k^ 



(4.5) 

Although the divergences can be deah with by dimensional regularization, this procedure 
is not well suited to numerical analysis (see however ref. [^). We will make our subtractions 
using an ultraviolet cutoff, Aa(t), constant in physical coordinates. This guarantees that the 
counterterms will be time independent. The renormalization then proceeds much in the 
same manner as in reference |T^; the quadratic divergences renormalize the mass and the 
logarithmic terms renormalize the quartic coupling and the coupling to the Ricci scalar. In 
addition, there is a quartic divergence which renormalizes the cosmological constant while 
the leading renormalizations of Newton's constant and the higher order curvature coupling 
are quadratic and logarithmic respectively. The renormalization conditions on the mass, 
coupling to the Ricci scalar and coupling constant are obtained from the requirement that 



the frequencies that appear in the mode equations are finite [13], i.e: 



A 



-ml + is-Rit) + '-fct>\t) + '-f{i>\t))B = -ml + ^nnt) + '-^cf\t) + ^(V^^(t))i^, 

(4.6) 

while the renormalizations of Newton's constant, the higher order curvature coupling, and 
the cosmological constant are given by the condition of finiteness of the semi-classical 
Einstein-Friedmann equation: 



A 



B 



A 



A 



R 



Finally we arrive at the following set of renormalizations ||T6|: 



SttNGr 

Or 
N 
Kr 
N 



SttNGb 
as 

N 



^R-l 



A' 



2 ln(A/K 



nil 



ln(A/fi:) 
IGvr^ 



IGvr^ 
A^ 



+ 



m|j ln(A//t) 



IGtt^ 



-m 



R 



:r 



+ A 



A2 



R 



, ^ . ln(A//€) 



167r2 

1\ ln(A//t 
" 6/ 167r2 



IGtt^ 



A 



mt)) 



R 



R 




e{k 



ka^{t) 



2k^ 



M\t) 



(4.7) 

(4.8) 
(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 
(4.14) 
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Here, k is the renormalization point. As expected, the logarithmic terms are consistent with 



the renormahzations found using dimensional regularization [pTiPq . Again, we set aji = 



and choose the renormalized cosmological constant such that the vacuum energy is zero in 
the true vacuum. We emphasize that while the regulator we have chosen does not respect 
the covariance of the theory, the renormalized energy momentum tensor defined in this way 
nevertheless retains the property of covariant conservation in the limit when the cutoff is 
taken to infinity. 

The logarithmic subtractions can be neglected because of the coupling A < 10~^^. Using 
the Planck scale as the cutoff and the inflaton mass m^j as a renormalization point, these 
terms are of order Aln[Mp;/m/j] < 10^^°, for m > 10^ GeV . An equivalent statement is 
that for these values of the coupling and inflaton masses, the Landau pole is well beyond the 
physical cutoff Mpi. Our relative error in the numerical analysis is of order 10~^, therefore our 
numerical study is insensitive to the logarithmic corrections. Though these corrections are 
fundamentally important, numerically they can be neglected. Therefore, in the numerical 
computations that follow, we will neglect logarithmic renormalization and subtract only 
quartic and quadratic divergences in the energy and pressure, and quadratic divergences in 
the fluctuation contribution. 



V. RENORMALIZED EQUATIONS OF MOTION FOR DYNAMICAL 

EVOLUTION 

It is convenient to introduce the following dimensionless quantities and definitions. 



H k 

T = niRt ; h = ; q = 



rriR 



87r2 ' 



(5.1) 



2m2 



02(t) ; ^S(r) 



2m% 



(5.2) 



Choosing ^i? = (minimal coupling) and the renormalization point k = \mR\ and setting 
a(0) = 1, the equations of motion become: 



ri{r) = 0, 



(5.3) 



l+^2^^S(r) 



= 0, 



/,(0) = [-/i(0)-za;,]/,(0). 



g2 + 1-7^2(0) + 



6m 



+ ^S(0) 



^S(O) 



for q^ > -1 +?7^(0) - 



7^(o) 



+ 9^0), 



for q' <-l+r]\0)--^ + gi:{0). 

brriR 



(5.4) 
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where 



S(r) 



q^dq 



1 , e(g-i) /,„ . .2 



a(r) 



The initial conditions for //(r) will be specified later. An important point to notice 
is that the equation of motion for the q = mode coincides with that of the zero mode 
(^.31). Furthermore, for rj^r oo) 7^ 0, a stationary (equilibrium) solution of the eq.(^.3|) is 
obtained when the sum rule ||T2|0,[16| 



5(S(oo) 







(5.5) 



is fulfilled. This sum rule is nothing but a consequence of Goldstone's theorem and is a 
result of the fact that the large N approximation satisfies the Ward identities associated 
with the 0{N) symmetry, since the term — 1 + 77^ + gTj is seen to be the effective mass of 
the modes transverse to the symmetry breaking direction, i.e. the Goldstone modes in the 
broken symmetry phase. 

In terms of the zero mode //(r) and the quantum mode function given by eq. (|5.4| ) we find 
that the Friedmann equation for the dynamics of the scale factor in dimensionless variables 
is given by 



hi 



(5.6) 



and the renormalized energy and pressure are given by: 



2 

9_ 
2 



/ q'dq \\f,\' - 5«(g, r) + ^ (|// - Q{q - 1) S'-'^q, r) 



(5.7) 



{P + ^)R{r) 



2Nm% 



A 



R 



+ 9 J q'dq \f,\' - S^'\q, r) + ^ {\f,\' - 0(g - 1) S^'\q, r 



(5.8) 



where the subtractions iS*-^-* and iS^^-* are given by the right hand sides of eqns. ([4.4|) and 
respectively. 

The renormalized energy and pressure are covariantly conserved: 

eRiT)+3h{r){p + e)niT) = 0. (5.9) 

In order to provide the full solution we now must provide the values of ri{0), rj{0), and 
hf). Assuming that the infiationary epoch is associated with a phase transition at the GUT 
scale, this requires that Nm^/Xn (10^^ Gev Y cind assuming the bound on the scalar 
self-coupling Xr 10^^^ — 10~^^ (this will be seen later to be a compatible requirement), we 
find that Hq ^ N^/'^ which we will take to be reasonably given by /iq ~ 1 — 10 (for example 
in popular GUT's ^ 20 depending on particular representations). 

We will begin by studying the case of most interest from the point of view of describing 
the phase transition: 77 (0) = and ?7(0) = 0, which are the initial conditions that led to 
puzzling questions. With these initial conditions, the evolution equation for the zero mode 
eq. (|5.3|) determines that //(r) = by symmetry. 
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A. Early time dynamics: 



Before engaging in the numerical study, it proves illuminating to obtain an estimate of 
the relevant time scales and an intuitive idea of the main features of the dynamics. Because 
the coupling is so weak {g ~ 10^^^ ^ 1) and after renormalization the contribution from 
the quantum fluctuations to the equations of motion is finite, we can neglect all the terms 
proportional to g in eqs.( |5.7|) and ( |5.4|) . 

For the case where we choose t]{t) =0 and the evolution equations for the mode functions 
are those for an inverted oscillator in De Sitter space-time, which have been studied by Guth 
and Pi POl 



One obtains the approximate solution 



h{t) ^ ho, 



An J,j I -: — e 





Bn J-u — e 



V 



19 1 
4 + ^' 



(5.10) 



where J±u{z) are Bessel functions, and Aq and Bg are determined by the initial conditions 
on the mode functions: 



A, 



1 



TT q 



^JuTqlhi^ sinz/TT 
1 TT g 



1 



'uJ^2ho sinz/TT 



(5.11) 
(5.12) 



After the physical wavevectors cross the horizon, i.e. when qe /h^ <C 1 we find that 
the mode functions factorize: 



Bn 



'2h 



Tl-u) 



0_\ g(z.-3/2)/ior_ 



(5.13) 



This result reveals a very important feature: because of the negative mass squared term 
in the matter Lagrangian leading to symmetry breaking (and u > 3/2), we see that all 
of the mode functions grow exponentially after horizon crossing (for positive mass squared 
u < 3/2, and they would decrease exponentially after horizon crossing). This exponential 
growth is a consequence of the spinodal instabilities which is a hallmark of the process of 
phase separation that occurs to complete the phase transition. We note, in addition that 
the time dependence is exactly given by that of the q = mode, i.e. the zero mode, which is 
a consequence of the redshifting of the wavevectors and the fact that after horizon crossing 
the contribution of the term q'^/a^{T) in the equations of motion become negligible. We 
clearly see that the quantum fluctuations grow exponentially and they will begin to be of 
the order of the tree level terms in the equations of motion when 5'S(r) 1. At large times 

S(r) ^ ^2(^^)^2g(2.-3)hor ^ 

with J^(/io) a finite constant that depends on the initial conditions and is found numerically 
to be of C(l) [see figure (0)]. 
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In terms of the initial dimensionful variables, the condition 5'S(r) ^ 1 translates to 
< ijj'^{x,t) >_R~ 2m\/\R, i.e. the quantum fluctuations sample the minima of the (renor- 
malized) tree level potential. We find that the time at which the contribution of the quantum 
fluctuations becomes of the same order as the tree level terms is estimated to be 



{2v - 3)ho 



In 



ghlJ="^{hQ 



-hn In 



1 



ghlJ^^ho) 



+ O{l/ho 



(5.14) 



At this time, the contribution of the quantum fluctuations makes the back reaction very im- 
portant and, as will be seen numerically, this translates into the fact that also determines 
the end of the De Sitter era and the end of inflation. The total number of e-folds during the 
stage of exponential expansion of the scale factor (constant ho) is given by 



2z/-3 



In 



g hi J="^{ho 



-hi In 
2 ° 



g hi J="^{ho 



Oil) 



(5.15) 



For large ho we see that the number of e-folds scales as hi as well as with the logarithm 
of the inverse coupling. These results ( |5.13| , p?T^ , p.l5D will be confirmed numerically below 
and will be of paramount importance for the interpretation of the main consequences of the 
dynamical evolution. 



B. r](0) 7^ 0: classical or quantum behavior? 

Above we have analyzed the situation when ri{0) = (or in dimensionful variables 
(f)(0) = 0). The typical analysis of inflaton dynamics in the literature involves the classical 
evolution of (j){t) with an initial condition in which 0(0) is very close to zero (i.e. the top 
of the potential hill) in the 'slow-roll' regime, for which (j) -C 3H(p. Thus, it is important 
to quantify the initial conditions on (f){t) for which the dynamics will be determined by 
the classical evolution of 0(t) and those for which the quantum fluctuations dominate the 
dynamics. We can provide a criterion to separate classical from quantum dynamics by 
analyzing the relevant time scales, estimated by neglecting non-linearities and backreaction 
effects. We consider the evolution of the zero mode in terms of dimensionless variables, and 
choose 77(0) 7^ and r)(0) = 0. (r)(0) 7^ simply corresponds to a shift in origin of time). 
We assume 7^(0)^ << 1 which is the relevant case where spinodal instabilities are important. 
We find 

ri{r) ^ r/(0) e^'^'i)'^"" . (5.16) 

The non-linearities will become important and eventually terminate inflation when 
//(r) ^ 1. This corresponds to a time scale given by 

^ In [1/^(0)] 



If Tc is much smaller than the spinodal time given by eq.( p.l4| ) then the c/asszca/ evolution 



of the zero mode will dominate the dynamics and the quantum fluctuations will not become 
very large, although they will still undergo spinodal growth. On the other hand, if Tc ^ 
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the quantum fluctuations will grow to be very large well before the zero mode reaches the 
non-linear regime. In this case the dynamics will be determined completely by the quantum 
fluctuations. Then the criterion for the classical or quantum dynamics is given by 

77(0) ^ ho =^ classical dynamics 

^(0) <^ \/g ho =^ quantum dynamics (5.18) 

or in terms of dimensionful variables 0(0) ^ Hq leads to classical dynamics and 0(0) ^ Hq 
leads to quantum dynamics. 

However, even when the classical evolution of the zero mode dominates the dynamics, 
the quantum fluctuations grow exponentially after horizon crossing unless the value of 0(t) 
is very close to the minimum of the tree level potential. In the large approximation the 
spinodal line, that is the values of 0(t) for which there are spinodal instabilities, reaches all 
the way to the minimum of the tree level potential as can be seen from the equations of 
motion for the mode functions. Therefore even in the classical case one must understand 
how to deal with quantum fluctuations that grow after horizon crossing. 



C. Numerics 

The time evolution is carried out by means of a fourth order Runge-Kutta routine with 
adaptive stepsizing while the momentum integrals are carried out using an 11-point Newton- 
Cotes integrator. The relative errors in both the differential equation and the integration are 
of order 10"*^. We flnd that the energy is covariantly conserved throughout the evolution to 
better than a part in a thousand. Figures (|l]-^ show 5'S(r) vs. r, /i(r) vs. r and ln(|/g(r)p) 
vs. T for several values of q with larger q's corresponding to successively lower curves. Figures 
(|,^ show p{T)/e{T) and the horizon size h'^ij) for g = 10~^^ ; 7^(0) = ; r]{0) = and we 
have chosen the representative value ho = 2.0. 

Figures (|l|) and (H) show clearly that when the contribution of the quantum fluctuations 
(y'S(r) becomes of order 1 inflation ends, and the time scale for 5'S(r) to reach 0{1) is very 
well described by the estimate (|5.14|) . From flgure 1 we see that this happens for r = ^ 90, 



leading to a number of e-folds A^e ~ 180 which is correctly estimated by (|5.14| , |5.15|) . 

Figure (^ shows clearly the factorization of the modes after they cross the horizon as 
described by eq. (|5.13|) . The slopes of all the curves after they become straight lines in flgure 
(0) is given exactly by (2z/ — 3), whereas the intercept depends on the initial condition on the 
mode function and the larger the value of q the smaller the intercept because the amplitude 
of the mode function is smaller initially. Although the intercept depends on the initial 
conditions on the long-wavelength modes, the slope is independent of the value of q and is 
the same as what would be obtained in the linear approximation for the square of the zero 
mode at times long enough that the decaying solution can be neglected but short enough 
that the effect of the non-linearities is very small. Notice from the flgure that when inflation 
ends and the non-linearities become important all of the modes effectively saturate. This 
is also what one would expect from the solution of the zero mode: exponential growth in 
early-intermediate times (neglecting the decaying solution), with a growth exponent given by 
(z/ — 3/2) and an asymptotic behavior of small oscillations around the equilibrium position, 
which for the zero mode is 77 = 1, but for the q ^ modes depends on the initial conditions. 
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All of the mode functions have this behavior once they cross the horizon. We have also 
studied the phases of the mode functions and we found that they freeze after horizon crossing 
in the sense that they become independent of time. This is natural since both the real and 
imaginary parts of /g(r) obey the same equation but with different boundary conditions. 
After the physical wavelength crosses the horizon, the dynamics is insensitive to the value 
of q for real and imaginary parts and the phases become independent of time. Again, this 
is a consequence of the factorization of the modes. 

The growth of the quantum fluctuations is sufficient to end inflation at a time given by 
in eq.( ^?T^ ). Furthermore figure shows that during the infiationary epoch p(r)/e(r) k, — 1 
and the end of infiation is rather sharp at with p{T)/e{T) oscillating between ±1 with zero 
average over the cycles, resulting in matter domination. Figure (|]) shows this feature very 
clearly; /i(r) is constant during the De Sitter epoch and becomes matter dominated after 
the end of infiation with h'^ij) ^ |(r — r^). There are small oscillations around this value 
because both p(r) and 6{t) oscillate. These oscillations are a result of small oscillations 
of the mode functions after they saturate, and are also a feature of the solution for a zero 
mode. 

All of these features hold for a variety of initial conditions. As an example, we show in 
figures (H) - d^) the plots corresponding to figs. (|1]) - (^) for the case of an initial Hubble 
parameter of /iq = 10. 



D. Zero Mode Assembly: 

This remarkable feature of factorization of the mode functions after horizon crossing can 
be elegantly summarized as 

fk{t)\k^nit)<^Ho = 9{q, ho)fo{r), (5.19) 

with kphit) = ke~^°^ being the physical momentum, g{q, ho) a complex constant, and /o(t) 
a real function of time that satisfies the mode equation with g = and real initial conditions 
which will be inferred later. Since the factor g{q, ho) depends solely on the initial conditions 
on the mode functions, it turns out that for two mode functions corresponding to momenta 
ki, ^2 that have crossed the horizon at times ti > t2, the ratio of the two mode functions at 
time t, (ts > t > ti > t2) is 



Then if we consider the contribution of these modes to the renorma/uec/ quantum fiuctuations 
a long time after the beginning of infiation (so as to neglect the decaying solutions), we find 
that 

where 'small' stands for the contribution of mode functions associated with momenta that 
have not yet crossed the horizon at time r, which give a perturbatively small (of order A) 
contribution. We find that several e-folds after the beginning of infiation but well before 
infiation ends, this factorization of superhorizon modes implies the following: 



16 



g I qHq |/^^(r) 
9 I q'dq \Bt) 



9 



a2(r) 



dq \mr) 



|CoP/J(r), 
IColVo'(r), 



(5.20) 
(5.21) 

(5.22) 



where we have neglected the weak time dependence arising from the perturbatively small 
contributions of the short- wavelength modes that have not yet crossed the horizon, and the 
integrals above are to be understood as the fully renormalized (subtracted), finite integrals. 
For ?7 = 0, we note that ( p.20[ ) and the fact that /o(t) obeys the equation of motion for the 
mode with q = leads at once to the conclusion that in this regime [(^^(r)]^ = |Co|/o(t) 
obeys the zero mode equation of motion 



|^ + 3/.A_i + (|Co|/o(.))^ 



\Co\foir) = 



(5.23) 



It is clear that several e-folds after the beginning of inflation, we can define an effective 
zero mode as 



Tjlffir) = gT,{r), or in dimensionful variables, 4>eff{t) = {i^'^{x,t)) ji 



(5.24) 



Although this identification seems natural, we emphasize that it is by no means a trivial or 
ad- hoc statement. There are several important features that allow an unambiguous identifi- 
cation: i) [{ip^{x,t))R] is a fully renormalized operator product and hence finite, ii) because 
of the factorization of the superhorizon modes that enter in the evaluation of [('^/'^(x, t))/j], 
^effi^) (|5.24|) obeys the equation of motion for the zero mode, iii) this identification is valid 
several e-folds after the beginning of inflation, after the transient decaying solutions have 
died away and the integral in {■ip'^{x,t)) is dominated by the modes with wavevector k that 
have crossed the horizon at t{k) -C t. Numerically we see that this identification holds 
throughout the dynamics except for a very few e-folds at the beginning of inflation. This 
factorization determines at once the initial conditions of the effective zero mode that can be 
extracted numerically: after the first few e-folds and long before the end of inflation we find 



<peff{t) = 0e//(O) e 



(5.25) 



where we parameterized 



ZTT 



to make contact with the literature. As is shown in fig. (pUj), we find numerically that 
J-'lHo/m) ^ 0{1) for a large range of 0.1 < Ho/m < 50 and that this quantity depends on 
the initial conditions of the long wavelength modes. 

Therefore, in summary, the effective composite zero mode obeys 



dr"^ dr 



3 

Veffir) = ; r]effir = 0) = (z/ - -) VeffiO) 



(5.26) 
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where ?7e//(0) = <^e//(0) is obtained numerically for a given by fitting the interme- 

diate time behavior of 5'S(r) with the growing zero mode solution. 

Furthermore, this analysis shows that in the case rj = 0, the renormalized energy and 
pressure in this regime in which the renormalized integrals are dominated by the superhori- 
zon modes are given by 

^«(^) - ^eff + i + ^eff) I , (5.27) 

iP + ^)n^^-^{nlff}, (5-28) 

where we have neglected the contribution proportional to l/a^(r) because it is effectively 
red-shifted away after just a few e-folds. We found numerically that this term is negligible 
after the interval of time necessary for the superhorizon modes to dominate the contribution 
to the integrals. Then the dynamics of the scale factor is given by 

h\r)=Ahl{lff^,, + \{-l+,!„y]. (5.29) 



We have numerically evolved the set of effective equations { f).2(: , 5.29| ) by extracting the 



initial condition for the effective zero mode from the intermediate time behavior of gT,{T). 
We found a remarkable agreement between the evolution of r^^jj and 5'S(r) and between the 
dynamics of the scale factor in terms of the evolution of ?7ejj(r), and the dynamics of the 
scale factor and quantum fluctuations within our numerical accuracy. Figures (pJ] ) and (|12|) 
show the evolution of ri1ff{T) and h{T) respectively from the classical evolution equations 
( p.2(j| ) and ( |5.29| ) using the initial condition rjefflO) extracted from the exponential fit of 



(7S(r) in the intermediate regime. These figures should be compared to Figure (|ID and 
We have also numerically compared p/e given solely by the dynamics of the effective zero 
mode and it is again numerically indistinguishable from that obtained with the full evolution 
of the mode functions. 

This is one of the main results of our work. In summary: the modes that become 
superhorizon sized and grow through the spinodal instabilities assemble themselves into an 
effective composite zero mode a few e-folds after the beginning of inflation. This effective 
zero mode drives the dynamics of the FRW scale factor, terminating inflation when the non- 
linearities become important. In terms of the underlying fluctuations, the spinodal growth of 
superhorizon modes gives a non-perturbatively large contribution to the energy momentum 
tensor that drives the dynamics of the scale factor. Inflation terminates when the mean 
square root fluctuation probes the equilibrium minima of the tree level potential. 

This phenomenon of zero mode assembly, i.e. the 'classicalization' of quantum mechan- 
ical fluctuations that grow after horizon crossing is very similar to the interpretation of 
'decoherence without decoherence' of Starobinsky and Polarski ||28|| . 

The extension of this analysis to the case for which 77(0) 7^ is straightforward. Since 



both ■r]{T) and ^ gT,{T) = |Co|/o(t) obey the equation for the zero mode, eq. (|5.3| ), it is clear 
that we can generalize our definition of the effective zero mode to be 



Veff[r) 



^r]^{T)+gi:{T) . (5.30) 
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which obeys the equation of motion of a classical zero mode: 



Veff{r)=0. (5.31) 

If this effective zero mode is to drive the FRW expansion, then the additional condition 

TiVo -2w/o/o + r/Vo' = (5.32) 

must also be satisfied. One can easily show that this relation is indeed satisfied if the mode 
functions factorize as in ( |5.19| ) and if the integrals ( |5.2CI| ) - ( |5.22| ) are dominated by the 



contributions of the superhorizon mode functions. This leads to the conclusion that the 
gravitational dynamics is given by eqns. ( [5.27] ) - ( |5.29| ) with f]eff{,T) defined by ( [5.301 ). 



We see that in all cases, the full large quantum dynamics in these models of inflation- 
ary phase transitions is well approximated by the equivalent dynamics of a homogeneous, 
classical scalar field with initial conditions on the effective field rieff{0) > y/ghoJ^{ho). We 
have verified these results numerically for the field and scale factor dynamics, finding that 
the effective classical dynamics reproduces the results of the full dynamics to within our 
numerical accuracy. We have also checked numerically that the estimate for the classical to 
quantum crossover given by eq. ( p. 18 ) is quantitatively correct. Thus in the classical case in 



which 77(0) ^ a/a ho we find that Veffi^) = vi'^)^ whereas in the opposite, quantum case 



Veffir) = ^gT.{T). 

This remarkable feature of zero mode assembly of long-wavelength, spinodally unstable 
modes is a consequence of the presence of the horizon. It also explains why, despite the fact 
that asymptotically the fluctuations sample the broken symmetry state, the equation of state 
is that of matter. Since the excitations in the broken symmetry state are massless Goldstone 
bosons one would expect radiation domination. However, the assembly phenomenon, i.e. 
the redshifting of the wave vectors, makes these modes behave exactly like zero momentum 
modes that give an equation of state of matter (upon averaging over the small oscillations 
around the minimum). 

Subhorizon modes at the end of inflation with q > ho e^°'^'' do not participate in the 
zero mode assembly. The behavior of such modes do depend on q after the end of inflation. 
Notice that these modes have extremely large comoving q since hoe''°^' > 10^6. As discussed 
in ref. |TH] such modes decrease with time after inflation as ~ l/a(r). 



VI. MAKING SENSE OF SMALL FLUCTUATIONS: 

Having recognized the effective classical variable that can be interpreted as the compo- 
nent of the field that drives the FRW background and rolls down the classical potential hill, 
we want to recognize unambiguously the small fluctuations. We have argued above that 
after horizon crossing, all of the mode functions evolve proportionally to the zero mode, and 
the question arises: which modes are assembled into the effective zero mode whose dynamics 
drives the evolution of the FRW scale factor and which modes are treated as perturbations? 
In principle every k mode provides some spatial inhomogeneity, and assembling these 
into an effective homogeneous zero mode seems in principle to do away with the very in- 
homogeneities that one wants to study. However, scales of cosmological importance today 
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first crossed the horizon during the last 60 or so e-folds of infiation. Recently Grishchuk 
29| has argued that the sensitivity of the measurements of AT/T probe inhomogeneities 



on scales ~ 500 times the size of the present horizon. Therefore scales that are larger than 
these and that have first crossed the horizon much earlier than the last 60 e-folds of infiation 
are unobservable today and can be treated as an effective homogeneous component, whereas 
the scales that can be probed experimentally via the CMB inhomogeneities today must be 
treated separately as part of the inhomogeneous perturbations of the CMB. 

Thus a consistent description of the dynamics in terms of an effective zero mode plus 
'small' quantum fluctuations can be given provided the following requirements are met: 
a) the total number of e-folds N^. ^ 60, b) all the modes that have crossed the horizon 
before the last 60-65 e-folds are assembled into an effective classical zero mode via 0e//(^) = 

[00 (^) + {"^"^{^yt)) r]^ , c) the modes that cross the horizon during the last 60-65 e-folds 
are accounted as 'small' perturbations. The reason for the requirement a) is that in the 
separation (p{x,t) = (peffit) + 6(j){x,t) one requires that 6(f){x,t)/(j)eff{t) -C 1. As argued 
above, after the modes cross the horizon, the ratio of amplitudes of the mode functions 
remains constant and given by e^'^'a)^^ with AA^ being the number of e-folds between the 
crossing of the smaller k and the crossing of the larger k. Then for S(j){x, t) to be much smaller 
than the effective zero mode, it must be that the Fourier components of 6(f) correspond to 
very large fc's at the beginning of inflation, so that the effective zero mode can grow for a 
long time before the components of S(f) begin to grow under the spinodal instabilities. In 
fact requirement a) is not very severe; in the flgures (1-5) we have taken ho = 2.0 which 
is a very moderate value and yet for A = 10~^^ the inflationary stage lasts for well over 
100 e-folds, and as argued above, the larger ho for flxed A, the longer is the inflationary 
stage. Therefore under this set of conditions, the classical dynamics of the effective zero 
mode (peffit) drives the FRW background, whereas the inhomogeneous fluctuations 6(j){x, t), 
which are made up of Fourier components with wavelengths that are much smaller than the 
horizon at the beginning of inflation and that cross the horizon during the last 60 e-folds, 
provide the inhomogeneities that seed density perturbations. 



VII. SCALAR AND TENSOR METRIC PERTURBATIONS: 



A. Scalar Perturbations: 



Having identifled the effective zero mode and the 'small perturbations', we are now in 
position to provide an estimate for the amplitude and spectrum of scalar metric pertur- 
bations. We use the clear formulation by Mukhanov, Feldman and Brandenberger 



m 



terms of gauge invariant variables. In particular we focus on the dynamics of the Bardeen 
potential which in longitudinal gauge is identifled with the Newtonian potential. The 
equation of motion for the Fourier components (in terms of comoving wavevectors) for this 
variable in terms of the effective zero mode is 1301 



H{t) - 



^ + 2fif(t)-i7(t)^' 



= 0. 



(7.1) 



We are interested in determining the dynamics of $fc for those wavevectors that cross 



20 



the horizon during the last 60 e-folds before the end of inflation. During the inflationary 
stage the numerical analysis yields to a very good approximation 



H{t)^Ho; 0e//(t) = 0e//(O) e^^^-i)^"*, 
where Hq is the value of the Hubble constant during inflation, leading to 



(1) 



f3 = u-l. 



(7.2) 



(7.3) 



The coefficients a^, hk are determined by the initial conditions. 

Since we are interested in the wavevectors that cross the horizon during the last 60 
e-folds, the consistency for the zero mode assembly and the interpretation of 'small pertur- 
bations' requires that there must be many e-folds before the last 60. We are then considering 
wavevectors that were deep inside the horizon at the onset of inflation. Mukhanov et. al. 
3Q| show that ^kif) is related to the canonical 'velocity field' that determines scalar per- 



turbations of the metric and which is quantized with Bunch-Davies initial conditions for the 
large fc-mode functions. The relation between $fc and v and the initial conditions on v lead 



at once to a determination of the coefficients and hk for k >> Hq 



mil 



^eff 



(0). 



TT 







(7.4) 



Thus we find that the amplitude of scalar metric perturbations after horizon crossing is 
jiven by 



(0) 



2^0 
k 



^{2u-3)Hot 



(7.5) 



The power spectrum per logarithmic k interval is given by |(5fc(t)p. The time dependence 
of |5fc(t)| displays the unstable growth associated with the spinodal instabilities of super- 
horizon modes and is a hallmark of the phase transition. This time dependence can be also 
understood from the constraint equation that relates the Bardeen potential to the gauge 
invariant field fluctuations [^, which in longitudinal gauge are identified with 6(j){x,t). 
The constraint equation and the evolution equations for the gauge invariant scalar field 
fluctuations are 



An 



(7.6) 



+ 3H- + — +M^ 



dt 



^eff 



^k + 2V'( 



= 



(7.7) 



Since the right hand side of ( [7.6|) is proportional to <peff/Mpi <C 1 during the inflationary 
epoch in this model, we can neglect the terms proportional to $fc and $fe on the left hand 
side of ( |7.7|) , in which case the equation for the gauge invariant scalar field fluctuation is 
the same as for the mode functions. In fact, since 6(j)f' is gauge invariant we can evaluate 
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it in the longitudinal gauge wherein it is identified with the mode functions /fc(t). Then 
absorbing a constant of integration in the initial conditions for the Bardeen variable, we find 



and using that oc ei'^-^/'^)^ot that after horizon crossing fk{t) oc e^''"^/^)-^"*, one 
obtains at once the time dependence of the Bardeen variable after horizon crossing. In 
particular the time dependence is found to be oc e^'^'^~^^^°*. It is then clear that the time 
dependence is a refiection of the spinodal (unstable) growth of the superhorizon field fiuc- 
tuations. 

To obtain the amplitude and spectrum of density perturbations at second horizon crossing 
we use the conservation law associated with the gauge invariant variable PD | 



3 1+ p/e 

which is valid after horizon crossing of the mode with wavevector k. Although this conserva- 
tion law is an exact statement of superhorizon mode solutions of eq.( [7.1|) , we have obtained 
solutions assuming that during the infiationary stage H is constant and have neglected the 
H term in Eq. (|7. 1| ) . Since during the infiationary stage, 

m = 0^,/t) oc H^, [^^J « Hi (7.10) 

and (p/cj) ~ ifo) the above approximation is justified. We then see that cfl^jit) oc e^'^'^^'^^^°^ 
which is the same time dependence as that of ^k{t)- Thus the term proportional to 1/(1 + 
p/e) in Eq. ( |7.9| ) is indeed constant in time after horizon crossing. On the other hand, 
the term that does not have this denominator evolves in time but is of order (1 + p/e) = 
—2H/3H^ -C 1 with respect to the constant term and therefore can be neglected. Thus, we 
confirm that the variable ^ is conserved up to the small term proportional to (1 +p/e)^k 
which is negligible during the infiationary stage. This small time dependence is consistent 
with the fact that we neglected the H term in the equation of motion for $fc(t). 

The validity of the conservation law has been recently studied and confirmed in different 
contexts |]32|,Q . Notice that we do not have to assume that $fc vanishes, which in fact does 
not occur. 

However, upon second horizon crossing it is straightforward to see that ^k{tf) ~ 0. The 
reason for this assertion can be seen as follows: eq.( [7.7|) shows that at long times, when the 
effective zero mode is oscillating around the minimum of the potential with a very small 
amplitude and when the time dependence of the fluctuations has saturated (see flgure 3), 
$fc will redshift as ~ l/a{t) and its derivative becomes extremely small. 

Using this conservation law, assuming matter domination at second horizon crossing. 



and ^k{tf) ~ we flnd 



'^^^'^)' = 5(.-|)^(i/o/-) ' ^'-''^ 
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where J-'{Ho/m) determines the initial amphtude of the effective zero mode ( |5.25|) . We can 
now read the power spectrum per logarithmic k interval 



Vsik) = \Sk\^ oc k 



leading to the index for scalar density perturbations 

n, = 1 - 2 ( z/ - ^ 



(7.12) 



(7.13) 



For Ho/m ^ 1, we can expand — 3/2 as a series in m? / Hq in eq. ( [7. 1 1| ) . Given that 
the comoving wavenumber of the mode which crosses the horizon n e-folds before the end of 
inflation is = Hqc^^"'^^ where Ne is given by ( |5.15| ), we arrive at the following expression 
for the amplitude of fluctuations on the scale corresponding to n in terms of the De Sitter 
Hubble constant and the couphng A: 



Mtf)\ 



_9Hl_ 

5y2m= 



(2e")"'/^^o VX 



2m^ /7 



3/fo V6 



-ln2 



+ 



m 



(7.14) 



Here, 7 is Euler's constant. Note the explicit dependence of the amplitude of density per- 
turbations on V^- For n ~ 60, the factor exp(nm^/3-ffQ) is (9(100) for Hq/itl = 2, while it 
is 0{1) for Hq/itl > 4. Notice that for Hq/iti large, the amplitude increases approximately 
as {Ho/m)^, which will place strong restrictions on A in such models. 

We remark that we have not included the small corrections to the dynamics of the 
effective zero mode and the scale factor arising from the non-linearities. We have found 
numerically that these nonlinearities are only signiflcant for the modes that cross about 60 
e-folds before the end of inflation for values of the Hubble parameter Hq/ttlji > 5. The effect 
of these non-linearities in the large limit is to slow somewhat the exponential growth of 
these modes, with the result of shifting the power spectrum closer to an exact Harrison- 
Zeldovich spectrum with Ug = 1. Since for Ho/m^i > 5 the power spectrum given by (|7.13 ) 
differs from one by at most a few percent, the effects of the non-linearities are expected to 
be observationally unimportant. The spectrum given by (|7. 1 1| ) is similar to that obtained 
in references P,pU| although the amplitude differs from that obtained there. In addition, we 
do not assume slow roll for which (i/ — |) <^ 1, although this would be the case if Ne 3> 60. 

We emphasize an important feature of the spectrum: it has more power at long wave- 
lengths because u — 3/2 > 0. This is recognized to be a consequence of the spinodal insta- 
bilities that result in the growth of long wavelength modes and therefore in more power for 
these modes. This seems to be a robust prediction of new inflationary scenarios in which the 
potential has negative second derivative in the region of fleld space that produces inflation. 

It is at this stage that we recognize the consistency of our approach for separating the 
composite effective zero mode from the small fluctuations. We have argued above that many 
more than 60 e-folds are required for consistency, and that the small fluctuations correspond 
to those modes that cross the horizon during the last 60 e-folds of the inflationary stage. 
For these modes H^/k = e-^ot*{k) ^j^g^-g t*{k) is the time since the beginning of inflation of 
horizon crossing of the mode with wavevector k. The scale that corresponds to the Hubble 
radius today Aq = 27r/fco is the flrst to cross during the last 60 or so e-folds before the 
end of inflation. Smaller scales today will correspond to k > ko at the onset of inflation 
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since they will cross the first horizon later and therefore will reenter earlier. The bound 
on \Sko\ oc AT/T < 10~^ on these scales provides a lower bound on the number of e-folds 
required for these type of models to be consistent: 

12 ln(z/-|) 
iVe>60 + ^ ^ 7.15 

where we have written the total number of e-folds as = Hq t*{kQ) + 60. This in turn can 
be translated into a bound on the coupling constant using the estimate given by eq. (^.15|) . 



The four year COBE DMR Sky Map |34| gives n ^ 1.2 ± 0.3 thus providing an upper 
bound on u 

0<z/-^<0.05 (7.16) 

corresponding to Hq > 2.6. We then find that these values of ho and A ~ 10"^^ — 10~^^ 
provide sufficient e-folds to satisfy the constraint for scalar density perturbations. 



B. Tensor Perturbations: 

The scalar field does not couple to the tensor (gravitational wave) modes directly, and 
the tensor perturbations are gauge invariant from the beginning. Their dynamical evolution 



is completely determined by the dynamics of the scale factor [p0| , |35| . Having established 
numerically that the inflationary epoch is characterized by H/Hq <^ 1 and that scales 
of cosmological interest cross the horizon during the stage in which this approximation is 
excellent, we can just borrow the known result for the power spectrum of gravitational waves 



produced during inflation extrapolated to the matter era |30,35 



(7,17) 



Thus the spectrum to this order is scale invariant (Harrison-Zeldovich) with an amplitude 
of the order m^/XMpi. Then, for values of m ^ 10^^ — 10^^ Gev and A ^ 10"^^ — 10"^'^ 
one finds that the amplitude is < 10^^° which is much smaller than the amplitude of scalar 
density perturbations. As usual the amplification of scalar perturbations is a consequence 
of the equation of state during the inflationary epoch. 

VIII. CONTACT WITH THE RECONSTRUCTION PROGRAM: 

The program of reconstruction of the inflationary potential seeks to establish a rela- 
tionship between features of the inflationary scalar potential and the spectrum of scalar and 
tensor perturbations. This program, in combination with measurements of scalar and tensor 
components either from refined measurements of temperature inhomogeneities of the CMB 
or through galaxy correlation functions will then offer a glimpse of the possible realization of 
the inflation Such a reconstruction program is based on the slow roll approximation 

and the spectral index of scalar and tensor perturbations are obtained in a perturbative 



expansion in the slow roll parameters |36-3 
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We can make contact with the reconstruction program by identifying above with our 0e// 
after the first few e-folds of inflation needed to assemble the effective zero mode from the 
quantum fluctuations. We have numerically established that for the weak scalar coupling 
required for the consistency of these models, the cosmologically interesting scales cross the 
horizon during the epoch in which H ^ Hq ; (p^jf ~ (i^ ~ 3/2) Hq 0e// ; V ~ ^r/ A ^ tiff- 
in this case we find 

v{<Peff) = -{^ - \) ; e(0e//) ^ 0{\) « r7(0e//). (8.3) 



With these identifications, and in the notation of |36,^ the reconstruction program 



predicts the index for scalar density perturbations given by 

3 



2 



0(A), (8.4) 



which coincides with the index for the power spectrum per logarithmic interval with 
\5k\ given by eq. (|7.11|) . We must note however that our treatment did not assume slow roll 



for which (z/ — |) <C 1. Our self-consistent, non-perturbative study of the dynamics plus the 
underlying requirements for the identification of a composite operator acting as an effective 
zero mode, validates the reconstruction program in weakly coupled new inflationary models. 



IX. DECOHERENCE: QUANTUM TO CLASSICAL TRANSITION DURING 

INFLATION 

An important aspect of cosmological perturbations is that they are of quantum origin 
but eventually they become classical as they are responsible for the small classical metric 
perturbations. This quantum to classical crossover is associated with a decoherence process 



and has received much attention 28 



Recent work on decoherence focused on the description of the evolution of the density 



matrix for a free scalar massless field that represents the "velocity field" pO[ associated 
with scalar density perturbations [^]. In this section we study the quantum to classical 
transition of superhorizon modes for the Bardeen variable by relating these to the field 
mode functions and analyzing the full time evolution of the density matrix of the matter 



field. This is accomplished with the identification given by equation (|7.8| ) which relates 
the mode functions of the Bardeen variable with those of the scalar field. This relation 
establishes that in the models under consideration the classicality of the Bardeen variable 
is determined by the classicality of the scalar field modes. 

In the situation under consideration, long-wavelength field modes become spinodally 
unstable and grow exponentially after horizon crossing. The factorization ( |5.13| ) results in 
the phases of these modes "freezing out" . This feature and the growth in amplitude entail 
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that these modes become classical. The relation ( [7.8|) in turn implies that these features 
also apply to the superhorizon modes of the Bardeen potential. 

Therefore we can address the quantum to classical transition of the Bardeen variable 
(gravitational potential) by analyzing the evolution of the density matrix for the matter 
field. 

To make contact with previous work ||28|j3^ we choose to study the evolution of the field 
density matrix in conformal time, although the same features will be displayed in comoving 
time. 

The metric in conformal time takes the form 

ds^ = C\r){dT^ -dx^). (9.1) 
Upon a conformal rescaling of the field 

l(x,t) = x(^,r)/c(r), (9.2) 

the action for a scalar field becomes, after an integration by parts and dropping a surface 
term 

S = j d'xdr[\{^f-\{Vxf-V{x)]. (9.3) 

with 

V(x) = C\T) V (^^) - c\r) g x^ (9-4) 

where TZ = 6C"{T)/C^{T) is the Ricci scalar, and primes stand for derivatives with respect 
to conformal time T (for more details see the appendix of ref. [16|). As we can see from 
eq.( p.3| ), the action takes the same form as in Minkowski space-time with a modified potential 

v(x). 

The conformal time Hamiltonian operator, which is the generator of translations in T, 
is given by 

Hr = Jd'x[\ul+l{Vxr + V{x)], (9.5) 

with being the canonical momentum conjugate to x, = x'- Separating the zero mode 
of the field x 

x(x,r) = xo(T)5,,i + ^(x,r), (9.6) 

and performing the large factorization on the fluctuations we find that the Hamiltonian 
becomes linear plus quadratic in the fluctuations, and similar to a Minkowski space-time 
Hamiltonian with a T dependent mass term given by 



^^ + (e-^)7^ + ^x^(T) + ^(f) 



(9.7) 
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We can now follow the steps and use the results of ref. for the conformal time 
evolution of the density matrix by setting a{t) = 1 in the proper equations of that reference 
and replacing the frequencies by 

u;l{T) = P + M\T) . (9.8) 

The expectation value in Eq. (|9.7|) and that of the energy momentum tensor are obtained 
in this T evolved density matrix. [As is clear, we obtain in this way the self-consistent 
dynamics in the curved cosmological background 

The time evolution of the kernels in the density matrix (see [0) is determined by the 
mode functions that obey 



+ k^ + M\r) 



Fk{r) = 0. 



The Wronskian of these mode functions 

W(F, F*) = 



F, F'* 



(9.9) 



(9.10) 



is a constant. It is natural to impose initial conditions such that at the initial T the density 
matrix describes a pure state which is the instantaneous ground state of the Hamiltonian 
at this initial time. This implies that the initial conditions of the mode functions F^iT) be 
chosen to be (see ITI 



Fk{To) 



l^kiXo) 

With such initial conditions, the Wronskian ( p.lO| ) takes the value 

>V(F, F*) = -2i . 



(9.11) 



(9.12) 



The Heisenberg field operators xi^^'F) and their canonical momenta n^(x, T) can now 
be expanded as: 



x{x, T) = 

n^(f,T) 



d^k 
(2^0^ 
d^k 
(27r)3/2 



d^,F,{r) + dl.F*{r) 



Ak-x 



^kFL{T) + dKF':{r) 



Ak-x 



(9.13) 
(9.14) 



with the time independent creation and annihilation operators and at obeying canonical 
commutation relations. Since the fluctuation fields in comoving and conformal time are 
related by a conformal rescaling given by eq. (|9.2|) it is straightforward to see that the mode 
functions in comoving time t are related to those in conformal time simply as 



fk{t) 



C{T)- 



(9.15) 



Therefore the initial conditions given in Eq. ( p.ll| ) on the conformal time mode functions 
and the choice a()f:o) = C(^) = 1 imply the initial conditions for the mode functions in 
comoving time given by Eq. ( p.24|) . 
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In the large or Hartree (also in the self-consistent one-loop) approximation, the den- 
sity matrix is Gaussian, and defined by a normalization factor, a complex covariance that 
determines the diagonal matrix elements, and a real covariance that determines the mixing 
in the Schrodinger representation as discussed in ref. |jl3| (and references therein). 

That is, the density matrix takes the form 

p[<f , ^, T] = n A4(r) exp {-^A,(T) %iT) . ff_^,{T) - ^A^T) %{T) ■ r/^.^T) 

k 

- B,{T) ff^,{T) . ff_^,{r) ■ - ^-^.^T))} ; , (9.16) 

f^^.iT) = - X.iX) 5{k) ; ~Wr)=Ur)-X,iT)5,,,5{k). 

T^%{T) is the Fourier transform of n^(T, x). This form of the density matrix is dictated by 
the hermiticity condition 

p[$,l',T]=p*[<l,$,r]; 

as a result of this, Bk{T) is real. The kernel BkiT) determines the amount of 'mixing' in 
the density matrix since if = 0, the density matrix corresponds to a pure state because it 
factorizes into a wave functional depending only on $(■) times its complex conjugate taken 
at $(■). This is the case under consideration, since the initial conditions correspond to a 
pure state and under time evolution the density matrix remains that of a pure state |]13| . 



In conformal time quantization and in the Schrodinger representation in which the field 
X is diagonal the conformal time evolution of the density matrix is via the conformal time 
Hamiltonian (|9.5| ). The evolution equations for the covariances is obtained from those given 
in ref. |13[ by setting a(t) = 1 and using the frequencies ujI{T) = k"^ + Ai"^ (T) . In particular, 
by setting the covariance of the diagonal elements (given by equation (2.20) in ||13|; see also 
equation (2.44) of |ll). 



(9.17) 



More explicitly 



Mk{T)=Uk{TQ) exp 



Aj.iT) = -Aiog|F,(T)| = -a{t)-a{t) |log|/.(t)| 



ARkiT) 



|F,(T)P a{tr\Mt)\ 



2 ' 



(9.18) 



where Arj. and Ajj^ are respectively the real and imaginary parts of Aj^ and we have used 
the value of the Wronskian ( |9.12| ) in evaluating ( |9.18| ). 

The coefficients Ak{T) and Mk{T) in the gaussian density matrix (|9.16|) are completely 
determined by the conformal mode functions Fk{T) (or alternatively the comoving time 
mode functions fkif))- 
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Let us study the time behavior of these coefficients. During inflation, a{t) ~ e'*"*, and 
the mode functions factorize after horizon crossing, and superhorizon modes grow in cosmic 
time as in Eq. ( p.l3| ): 



where the coefficient can be read from eq. (|5.13| ) 



We emphasize that this is a result of the full evolution as displayed from the numerical so- 
lution in fig. (^. These mode functions encode all of the self-consistent and non-perturbative 
features of the dynamics. This should be contrasted with other studies in which typically 
free field modes in a background metric are used. 

Inserting this expression in eqs. (|9.18|) yields. 



AjkiT) e'^»* - ^) + , 

Since z/ — ^ > 1, we see that the imaginary part of the covariance Ajk{T) grows very fast. 
Hence, the off-diagonal elements of p[$, $, T] oscillate wildly after a few e-folds of inflation. 
In particular their contribution to expectation values of operators will be washed out. That 
is, we quickly reach a classical regime where only the diagonal part of the density matrix is 
relevant: 

p[$, $, T] = nAr.(T) exp {-A«.(T) r/g(T) r/_,,(T)} . (9.19) 

it 

The real part of the covariance Ajik{T) (as well as any non-zero mixing kernel Bk{T) 
||13|| ) decreases as e~(^'^~^)'*o*. Therefore, characteristic fleld conflgurations r]^ are very large 
(of order e^'^~5)'*o*). Therefore conflgurations with fleld amplitudes up to 0{e^'^~^^'^°^) will 
have a substantial probability of occurring and being represented in the density matrix. 

Notice that x ~ e*-'^~2'"*°* corresponds to fleld conflgurations $ with amplitudes of order 
g(!^-2)'^o« i^ggg ( p.2|) ]. It is the fact that z/ — | > which in this situation is responsible 
for the "classicalization" , which is seen to be a consequence of the spinodal growth of long- 
wavelength fluctuations. 

The equal-time fleld correlator is given by 

(X(f, T) x(x , T)) = / ^ \F,{T)\' e^-^-''^ , 

= <tf jl^,\m?e^'-''-''' ■ (9-20) 

and is seen to be dominated by the superhorizon mode functions and to grow as ei'^^-^)^ot^ 
whereas the fleld commutators remain flxed showing the emergence of a classical behavior. 
As a result we obtain 

(X(x,r) x{S\T)) oc a\t) (f),ff{t) (j),ff{t') G{\x~x'\) + small (9.21) 
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where G{\x — x'\) falls off exponentially for distances larger than the horizon [|I4| and "small" 
refers to terms that are smaller in magnitude. This factorization of the correlation functions 
is another indication of classicality. 

Therefore, it is possible to describe the physics by using classical field theory. More 
precisely, one can use a classical statistical (or stochastic) field theory described by the 
functional probability distribution (|9.19|) . 

These results generalize the decoherence treatment given in ref. for a free massless 
field in pure quantum states to the case of interacting fields with broken symmetry. Note that 
the formal decoherence or classicalization in the density matrix appears after the modes with 
wave vector k become superhorizon sized i.e. when the factorization of the mode functions 
becomes effective. 



X. CONCLUSIONS: 

It can be argued that the inflationary paradigm as currently understood is one of the 
greatest applications of quantum field theory. The imprint of quantum mechanics is every- 
where, from the dynamics of the inflaton, to the generation of metric perturbations, through 
to the reheating of the universe. It is clear then that we need to understand the quantum 
mechanics of inflation in as deep a manner as possible so as to be able to understand what 
we are actually testing via the CMBR temperature anisotropics, say. 

What we have found in our work is that the quantum mechanics of inflation is extremely 
subtle. We now understand that it involves both non-equilibrium as well as non-perturbative 
dynamics and that what you start from may not be what you wind up with at the end! 

In particular, we see now that the correct interpretation of the non-perturbative growth 
of quantum fluctuations via spinodal decomposition is that the background zero mode must 
be redefined through the process of zero mode reassembly that we have discovered. When 
this is done (and only when!) we can interpret inflation in terms of the usual slow- roll 
approach with the now small quantum fluctuations around the redefined zero mode driving 
the generation of metric perturbations. 

We have studied the non-equilibrium dynamics of a 'new inflation' scenario in a self- 
consistent, non-perturbative framework based on a large N expansion, including the dy- 
namics of the scale factor and backreaction of quantum fluctuations. Quantum fluctuations 
associated with superhorizon modes grow exponentially as a result of the spinodal insta- 
bilities and contribute to the energy momentum tensor in such a way as to end inflation 
consistently. 

Analytical and numerical estimates have been provided that establish the regime of 
validity of the classical approach. We find that these superhorizon modes re-assemble into 
an effective zero mode and unambiguously identify the composite field that can be used 
as an effective expectation value of the inflaton field whose classical dynamics drives the 
evolution of the scale factor. This identification also provides the initial condition for this 
effective zero mode. 

A consistent criterion is provided to extract "small" fluctuations that will contribute to 
cosmological perturbations from "large" non-perturbative spinodal fluctuations. This is an 
important ingredient for a consistent calculation and interpretation of cosmological pertur- 
bations. This criterion requires that the model must provide many more than 60 e-folds 
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to identify the 'small perturbations' that give rise to scalar metric (curvature) perturba- 
tions. We then use this criterion combined with the gauge invariant approach to obtain the 
dynamics of the Bardeen variable and the spectrum for scalar perturbations. 

We find that during the inflationary epoch, superhorizon modes of the Bardeen poten- 
tial grow exponentially in time reflecting the spinodal instabilities. These long-wavelength 
instabilities are manifest in the spectrum of scalar density perturbations and result in an 
index that is less than one, i.e. a "red" power spectrum, providing more power at long 
wavelength. We argue that this 'red' spectrum is a robust feature of potentials that lead 
to spinodal instabilities in the region in field space associated with infiation and can be 
interpreted as an "imprint" of the phase transition on the cosmological background. Tensor 
perturbations on the other hand, are not modified by these features, they have much smaller 
amplitude and a Harrison-Zeldovich spectrum. 

We made contact with the reconstruction program and validated the results for these 
type of models based on the slow-roll assumption, despite the fact that our study does not 
involve such an approximation and is non-perturbative. 

Finally we have studied the quantum to classical crossover and decoherence of quantum 
fiuctuations by studying the full evolution of the density matrix, thus making contact with 
the concept of "decoherence without decoherence" which is generalized to the interact- 
ing case. In the case under consideration decoherence and classicalization is a consequence 
of spinodal growth of superhorizon modes and the presence of a horizon. The phases of 
the mode functions "freeze out" and the amplitudes of the superhorizon modes grow ex- 
ponentially during the infiationary stage, again as a result of long-wavelength instabilities. 
As a result field configurations with large amplitudes have non-vanishing probabilities to be 
represented in the dynamical density matrix. In the situation considered, the quantum to 
classical crossover of cosmological perturbations is directly related to the "classicalization" 
of superhorizon matter field modes that grow exponentially upon horizon crossing during 
infiation. The diagonal elements of the density matrix in the Schroedinger representation 
can be interpreted as a classical distribution function, whereas the off-diagonal elements are 
strongly suppressed during inflation. 
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FIGURES 
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FIG. 1. c/S vs. r, for ri{0) = 0, 77(0) = 0,g = lO'^"^, ho = 2.0. 
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FIG. 2. i7(T) vs. T, for 77(0) = 0,?)(0) = 0,^ = 10~i4,/io = 2.0. 
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Figure 3 
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FIG. 3. ln(j/y(T)|2) vs. r, for r?(0) = 0,?7(0) = 0,g = 10''^^, ho = 2.0 for q = 0.0,5,10,15,20 
with smaller q corresponding to larger values of ln(|/g(r)p). 
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FIG. 4. p/e vs. T, for r]{0) = 0,77(0) = 0,g = 10~^^,ho = 2.0. 
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FIG. 5. 1//i(t) vs. r, for ri{0) = 0,?7(0) = 0,^ = W-^'^,ho = 2.0 
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FIG. 6. vs. T, for r/(0) = 0, 77(0) =0,g = IQ-^*, ^0 = 10.0. 
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Figure 7 
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FIG. 7. H{t) vs. r, for 7?(0) = 0,?7(0) = 0,^ = W^^^ho = 10.0. 
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FIG. 8. ln(|/g(r)p) vs. r, for r/(0) = 0,?7(0) = 0,g = 10"^^, /iq = 10.0 for q = 0.0,5,10,15,20 
with smaller q corresponding to larger values of ln(|/q(T)p). 
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Figure 9 
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FIG. 9. vie vs. r, for r/(0) = 0,7?(0) = 0,c/ = \^-^'^ M = 10.0 
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FIG. 10. T{H/m) vs. H, where J^{H/m) is defined by the relation (0) = {H/27r)T{H/m) 
[see eqs. (1^) and (HI)). 
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FIG. 11. vlffir) vs. r, for r/e//(0) = 3.94 x 10-^7?e//(0) = 0.31777e//(0), 5 = IQ-^^ho = 2.0. 
The initial conditions were obtained by fitting the intermediate time regime of 5S(t) in fig. (|l|). 
r/e//(r) is the solution of eq.( 5.26| ) with these initial conditions. 
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FIG. 12. /i(r) vs. T, obtained from the solution of eqns. ( ^.26| ) and ( ^.29 ) with the conditions 

of fig. (ITTD 
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